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Q ' We propose a scheme for generation of continuous-wave THz radiation. The 



scheme requires a medium where three discrete states in a A configuration can 
be selected, with the THz-frequency transition between the two lower metastable 
states. We consider the propagation of three-frequency continuous-wave electromag- 



> 

' netic (e.m.) radiation through a A medium. Under resonant excitation, the medium 



absorption can be strongly reduced due to quantum interference of transitions, while 

OS ; 

the nonlinear susceptibility is enhanced. This leads to very efficient energy transfer 
between the e.m. waves providing a possibility for THz generation. We demonstrate 
that the photon conversion efficiency is approaching unity in this technique. 



OCIS codes: 190.2620, 190.4410, 270.1670, 020.1670 

Generation of coherent terahertz electromagnetic radiation is a subject of much cur- 
rent research. This radiation might have interesting potential applications for electronics, 
chemical analysis of materials, local radars, tomography, environment monitoring^, and also 
for quantum optics applications and for frequency standards^!. The THz range is in the 
gap between the radio-frequency and the visible ranges. Therefore, neither conventional 
electronic (microwave) methods nor the photonic ones can be directly applied to generate 
coherent THz radiation. Pulsed THz fields are produced since early seventies and success- 
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fully used in different applications. However, to the best of our knowledge, a reliable source 
of continuous-wave (c.w.) THz radiation is not available up to now. 

Here we propose a c.w. generation scheme based on the effect of quantum interference in 
multilevel quantum systems (atoms, molecules, dopants in solids) induced by applied elec- 
tromagnetic radiation. In an optically dense media, the interference may be destructive for 
linear susceptibility (absorption and refraction), so that an otherwise opaque medium be- 
comes transparent. This process is termed electromagnetically induced transparency (EIT)@. 




The cancellation of absorption and refraction can also be explained as being due to the 
preparation of atoms in a coherent superposition ("dark" state) which is immune to applied 
radiation. When all the atomic population is perfectly trapped in this superposition, the 
medium does not "see" the radiation. It turns out, however, that if the dark state is slightly 
disturbed the linear susceptibility remains still very small, while the nonlinear suscepti- 
bility is resonantly enhanced by constructive interference!. Therefore, nonlinear frequency 
conversion/generation processes are very efficient in such media. Several successful experi- 
ments have been performed demonstrating efficient generation of XUV radiation with atomic 
hydrogen!, red to blue frequency conversion with molecular sodiumi, enhanced four-wave 
mixing with doped crystals^. Harris with coworkers have reached blue to UV§ and UV to 
VUV0 conversion in atomic Pb vapor with almost unity photon-conversion efficiency. The 
aim of the present paper is to extend these ideas to the case of THz generation, treating an 
interaction of the e.m. radiation with atoms as well as the propagation of radiation through 
the medium in exact manner. 

We suggest that the THz radiation can be produced by the EIT-assisted frequency con- 
version in media where three discrete states may be selected, two of them being metastable 
ones. An example of such a scheme is a closed three-level A system (Fig. 1), where |1) — |3) 
and |2) — 13) are the dipole-allowed optical transitions, and the frequency of a magnetic-dipole 
transition |1) — 12) is in the THz range. Such systems can be found in vapors of alkaline-earth 
atoms with the corresponding states |1) = 3 P or 3 P2, 1 2) = 3 -Pi and |3) = 3 Si. Frequencies 
of the transitions 3 P - 3 Pi and 3 P 2 - 3 Pi are 0.6 and 1.2 THz for 24 Mg, 1.5 and 3.2 THz 
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for 40 Ca, etc. A similar A system can be excited, for instance, in samarium atoms where the 
transition between two substates 4/ 6 6s 2 7 F and 7 Fi of the ground state fine structure has 
a frequency of 8.78 THz. One may also find appropriate schemes in other atoms, where the 
THz transition is usually between the two fine structure components, and probably in some 
molecules and doped crystals as well. 

Let us assume that the A atom interacts with a bichromatic optical field 

E(z,t) = ^2 e 3m E 3m (z,t)-exp[-i(uj 3m t - k 3m z + ip 3m (z,t))] + c.c. (1) 

m=l,2 ^ 

having the frequencies u 3m resonant or near-resonant with transitions \m) — |3) , m — 1, 2, 
and with the radiation of THz frequency lot 

H(z, t) = e T H(z, t)^- exp [—i (uj T t - k T z + ipr(z, t))] + c.c, (2) 

resonant with transition |1) — |2). Here, E 3m and H are the amplitudes, e 3m and are the 
unit polarization vectors, <p 3m and (p? are the phases, and k 3m = Lj 3m /c and kx = u)t/c are 
the wavenumbers of electric and magnetic field, respectively. Both the amplitudes E 3m , H 
and the phases (f 3m , <Pt are regarded as slowly varying functions of time and coordinate. 

The propagation of e.m. waves along the z axis in the medium is governed by the Maxwell 
equations. In the slowly varying amplitude and phase approximation, these equations can 
be transformed to the followine; 

dE 3m ldE 3m A A^d 3m uj 3m 

~dz~ + c-dT = ~ N ^^- Im(a3m) ' (3a) 

+ TW = ~ n ^^e\^ r ^ (3b) 

dH I OH 4vr /juj Tt . 

■^ + -cW = - N ^7 LlmM - (3c) 

-dz- + c^r = - N —H ReM > (3d) 

where d 3m = \e 3m d 3m \ , \i = |e^/2| , d 3m and ft are the matrix elements of the electric-dipole 
moment operator d and the magnetic-dipole moment operator /t, respectively, in the basis 
of bare atomic states \n), n = 1,2,3: d 3m = (3|d|m), jl = (l\p,\2); N is the density of the 
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active atoms. The medium polarization components (the right-hand side of Eqs. (|3])) are 
determined by the density matrix elements a ns averaged over the atomic velocities with the 
distribution w (v z ) where v z is the z-projection of the atom velocity: 



+00 



' ns 



dv z w (v z ) a ns , (4) 

—00 
with 

Ons = Pns (V z ) exp [l (u ns t - k ns Z + Xns)} , (5) 

where p ns = (n\p\s), p is the atomic density matrix; and the phase Xns is the sum of the 
e.m. field phase (p ns and the phase $ ns of the atomic dipole moment d 3m = Ic^le 1 ^ 3 ™, 

fl = \f2\ e 1 ^ 12 : X3m = <p3m + ^3m, Xl2 = + $12- 

The A system represents a typical interaction scheme where the interference of excitation 
channels may lead to creation of a dark state and preparation of atoms in this stateS If 
the e.m. field between the states |1) and \2) is not applied then the following superposition 
of the metastable states 

\NC) = 932 /931 |1> - exp ( X32 - X3 i) , l == |2> , (6) 

V 1 +^32/^31 V 1 +^32/531 

gets completely decoupled from the field when the two-photon resonance condition is satis- 
fied: 

A 32 - A 3 i = 0. (7) 

In Eqs. (H) and (0), g 3m = d 3m E 3m /2h are the Rabi frequencies of the optical fields, A 3m = 
w 3m — (£3 — Em) /% are the laser frequency detunings from transitions \m) — |3) , ( m = 1, 2), 
S n is the eigenenergy of the atomic state |n). The superposition \NC) is stable, and it is fed 
by spontaneous emission from the level |3) so that all the atomic population is trapped in 
\NC) after some optical pumping time. Therefore, the light ceases to interact with atoms 
and propagates through the medium without absorption and refraction. We note that this 
occurs for any relation between the laser intensities and phases. 



At the same time, the preparation of atoms in the superposition of metastable states |1) 
and |2) means that the coherence <Ti 2 is induced. Corresponding to the Eqs. (| c,d), this 
should lead to a generation of the field with frequency (o>3i — W32). Such an EIT-assisted 
generation of a microwave radiation has recently been observed in atomic Cs vapor^l. In 
general, the presence of a third e.m. field destroys the dark state, so that the optical waves 
start to interact with the medium. However, the degree of the destruction is very small 
if the generated field is weak. In the present case, the THz generation corresponds to the 
difference-frequency generation scheme^. If we suppose the photon conversion efficiency 
to be unity, that is for each photon of frequency U31 one THz-photon (and one photon of 
frequency U32) is generated, then the maximum energy (or power) conversion efficiency is 
given by 77 = huix/huj^i <C 1. For example, for the A system (3 3 Pi — 3 3 P2 — 4 3 Si) in 24 Mg, 
we have u T = 1.22 • 10 12 Hz, u 31 = 5.80 • 10 14 Hz, and r] = 0.0021. Therefore, the THz 
intensity is always much smaller than the optical one, so that the dark state disturbance is 
always very small. Thus, the total e.m. energy dissipation in the medium is negligible, and 
one can hope to get a very efficient generation in such a nearly dark ("grey") medium. 

In this paper we concentrate on a regime of the radiation propagation in a continuous 
wave limit, that is, we shall assume that the characteristic time of a change in the field 
amplitudes and phases, and the interaction time of atoms with the radiation are much longer 
than the characteristic time of a change in the internal state of an atom. Then, the time 
derivatives in Eqs. ([3]) can be dropped, and the steady-state values of the density matrix 
elements a ns can be used in the right-hand side of Eqs. (|3|). The system considered here is 
an example of the so-called closed-loop system, where the steady state can be established 
only when the multi-photon resonance condition is satisfied0: 

W31 - ^32 - uj t — 0. (8) 

In what follows we will always assume the condition (|j) to be fulfilled. 

The density matrix equations for the A atom can be found in Refs.00. These equations 
can be solved analytically under some simplifying conditions. We suppose that spontaneous 



relaxation rates for the channels |3) — > |1) and |3) — > |2) are equal: 731 = 732 = 7. With 
the condition cj 31 ss W32 (since (cj 31 — CJ32) /cj 31 = ujt/uj^ <C 1) this would also imply 
|d 3 i| ss |d 32 |. This simplification allows us to get analytical results which demonstrate all 
the basic features of the process. Numerical calculations with real atomic parameters (Figs. 
2-4) confirm the analytical results not only qualitatively but even quantitatively. This 
suggests that the performance of the proposed scheme is not very sensitive to exact values 
of the relaxation rates and dipole moments on optical transitions. Another assumption used 
is quite reasonable. Since the THz radiation intensity is very small, we assume that its Rabi 
frequency is small: Qt = pH/2h <C 7. 

Solution of the steady-state density matrix equations for v z = 0, T = (r is the relaxation 
rate of the coherence between states |1) and |2)), equal spontaneous relaxation rates: 731 = 
7 32 = 7, and zero detunings A 31 = A 32 = 0, gives to the first order in (g T / 7): 

Im ( 0-31 ) = sin $, (9a) 

9o 

Im(<732) = -^sin$, (9b) 
9o 

Im(<7 2 i) = -^f^sin$, (9c) 
9o 



Re(a 31 ) = - ^2(9h -9l) ^ (1Qa) 
9o 

Re(a 32 ) = 9T9 ' 1 {9h A - 9ll) cos $, (10b) 
9o 

Re(a 2 i) = C os $, (10c) 

9o 

The populations of the metastable states are 

gh 2gT93i932i . A 

pn = — 4 sin $, 

9o 9o 

#31 . 2 9T93l932l . ^ 

P22 = —o H 4 sin $, 

#0 9o 

where g$ = g\ x + g\ 2 i and the relative phase <3> is determined as 

$ = (X31 - X32) - Xl2- (11) 



The excited state population p 33 , which is responsible for irreversible dissipation of the e.m. 
energy by the medium, is of the second order in (g? / 7). This indicates that atoms are 
really in the grey state. One sees from Eqs. (|9]) and (|T0D, that the medium is absolutely 
transparent and not refractive for 

$ = 7m, n = 0,1,2,... (12) 

and 

031 = 932- (13) 



These are exactly the conditions for the dark state in closed A systemli3'EIrE 

For arbitrary optical field amplitudes and phases, however, the refraction and absorption 
(or amplification) of individual frequency components may be substantial. The change of 
the fields can be calculated analytically by the method developed in paper by Armstrong 
et. a/S. We insert the density matrix elements @ and (|I0| ) in the Maxwell equations @. 
The amplitude equations are read then: 
dE 31 nN d 31 d 32 /j 



dz h 2 c go 
dE 32 _ ttN d 31 d 3 2 fi 
dz h 2 c go 
dH ttN d 3 id 3 2H 



-lu 31 E 32 H sin $, (14a) 
u 32 E 31 H sin&, (14b) 
uutE 3 iE 3 2 sin (14c) 



dz h 2 c gl 

The intensities associated with each of these waves are given by I 3m = (c/87r) E\ m and Jy = 
(c/87r)H 2 . The set of equations ([14]) shows that the total power flow (proportional to the 
total intensity / = I 3 i + I 32 +It) is conserved, as expected for propagation through a lossless 
(due to EIT) medium: 

dl dl 31 dl 32 dI T N d 31 d 32f i 

Tz = llz~ + lLz~ + llz~ = i^ + uj T -uj 31 )E 31 E 32 Hsm^ = 0, 

where the last equality follows from the multiphoton resonance condition (||). Obviously, 
the energy losses (proportional to the population of the excited state p 33 ) appear only in the 
second order in (g? / 7). 



The set of equations (|14D also implies that 

A ( Isl ) = _A ( 732 A = _A ( It \ (15) 
dz \huj 3 iJ dz \hLU32J dz \%uj T ) ' 

which are the well-known Manley-Rowe relations^. These relations tell us that the rate at 
which photons at frequency uj? are created is equal to the rate at which photons at frequency 
u 32 are created and is equal to the rate at which photons at frequency u 32 are destroyed. 

In order to solve the Maxwell equations, we introduce the following dimensionless vari- 
ables: field amplitudes 

u m = — I E 3m , (m = 1, 2), 

V c uo J 

/8tt ujtY 1 ' 2 

\ c u J 

where ujq = (^31 + ^32 + ^r) /2, and a dimensionless optical length ( = nz with the coeffi- 
cient 

4nNd 31 d 32 fi /8tt / N 1 / 2 

K = = ^7^31^32- 

C g$ \ C UJqJ 

As we have discussed above, the THz intensity is always much smaller than the optical 
one: It <C /31 + / 3 2- Therefore, we can safely assume that the total optical intensity 
J 3 i + 1 32 ~ I is conserved. This follows also from the first of the Manley-Rowe relations flT5|) 
taking into account the relation u 3 i ~ u 32 ~ ujq (one can infer that an error in making such 
an approximation is of the order of 77 = (co> 31 — u 32 ) /uj 3 \ = uj t /oj 3 i < 1). The conservation 
of the optical intensity may be written in variables follows: 

«f(0 + «i(o = 1- (is) 

Other "constants of motion" following from the Manley-Rowe relations ([TJ) are: 

«i(0 + u 2 T (0 = B, (17) 
u|(0-^(C) = C (18) 

where constants B and C are determined from the boundary conditions at the ( = 0. 



Taking into account the assumed above closed values of the optical dipole moments: 
d 3 i « d 32 , one concludes that the quantity g\ is also approximately constant over all the 
optical path, and it is approximately equal to {2i[d\ x /fi 2 cj I. All this allows one to write 
the amplitude and the phase Maxwell equations in variables u m , $ as follows: 
dui 

— — = — u<iUt sin $, (19a) 
at, 

du 2 . _ , 1rin , 

— — = Ui%sm <P, (19b) 
ac, 

MiM 2 sm?>, (19c) 



d( 

rf$ My (ul — uf) — u\u\ 



cos$, (19d) 



d( U\U2Ut 

where we have assumed that the atomic dipole phases "& ns do not change along the propa- 
gation path. The optical length can now be re-written as: 

> irN /8tt \-^ 2 ; . , 

C = -a-/i — I y^TLJoZ. (20) 



he \ c 

The constants (|16D , ( fL7|) reduce the problem to solve of the set consisting of the first and 
the last of Eqs. fll9|) . The right-hand side of Eq. (]19| d) can be transformed to the following 
form by use of the first three Eqs. (HIT) : 



U% (Uo — uf) — u\u% COS$ 1 d(uiU2UT) 

cos <P — 



UiU 2 Ut sin $ UiU 2 Ut d( 

Therefore, the equation for the phase can be re-written as 

d$ cos$ d , 

-77 = — -rr^r In [uiU 2 Ut)) , 

dQ sm <P aC 

which can immediately be integrated to give the fourth constant of motion: 

uim 2 mtcos$ = II. (21) 

The value of the constant II can be determined from the known values of «j and $ at 
the entrance to the medium, £ = 0. The expression (|2"ID is used then to express sin $ = 



\J\ — cos 2 <& in terms of the conserved quantity II, and substitute it, together with the 
constants given in Egs-flHf), (0), into the Eq. ([19] a). This gives: 



,2 



C = ± I [ <Ml (22 ) 

2 ulL) y/uUl-ul)(B-ul)-IP 

The integral in Eq. fl22|) can be reduced to the elliptic integral. In this paper we consider 

the case of the THz radiation generation. This situation corresponds to wt(C = 0) = 0, B = 

uf = «i(C = 0), II = 0, and to the constant cos $: 

cos$(C) = 0. (23) 

Introduction of a quantity y = ui/u w leads to the elliptic integral of the first kind§ 



in 



standard form: 
v 



c = ±- * 



,„ V(l-» 2 )(l-«?o» 2 )' 
The solution can be written then through sinus Jacobi elliptic functions: 

uKO = u 2 w sn 2 [{( + ( );u w }, (24a) 
ul{C) = 1 - u 2 10 sn 2 [(C + Co) ; uio] , (24b) 
4(C) = «? (l - sn 2 [(( + Co) ; «io]) • (24c) 

The solution indicates that as the optical length increases, energy is periodically transferred 
between the e.m. waves. The initial condition ut{( = 0) = requires that the constant Co 
is equal to ±K(w 10 ), the complete elliptic integral 

dy 



K(u 



10 J 



{ y/(l-V 2 )(l-v.W)' 
which is a quarter-period of the function sn. Thus, the period of intensity oscillations is equal 

to 2K(«io) and the maximum possible power transferred to the THz radiation: u% = u\ 

occurs in a length 

Cmax = K(Uio). (25) 

The value of the function K(«io) is close to 7r/2 at Uio 1 and it is increased with u\q up 
to infinity at the limit u\q = 10. For u^q = 1 — u 2 Q < 1 a half period can be approximated 
by 

10 



(26) 



In terms of the real intensity, the Eq. ([24] c) for the THz radiation can be written: 



MO =MC = 0)— (l ~cd 2 [(;u l0 \). 



(27) 



We observe that at the length £ max = K(m 10 ) the optical field is converted into the THz 
radiation with unity photon efficiency. 

The obtained results are confirmed by numerical calculations of the Maxwell equations 
(§) with full density matrix calculations without any approximation. Figure 2 demonstrates 
the evolution of the e.m. wave intensities and the relative phase $ with the length for real 
parameters of the 3 3 Pi — 3 3 P2 — 4 3 Si system in 24 Mg atomic vapor. The optical length in 
this and all following figures is plotted in terms of a single-atom absorption cross-section for 
the optical field r = (3irc 2 /2uj 2 i ^J Nz, which has an advantage of being independent of total 
intensity /. The length ( introduced in Eq. ([20D is connected to r by 



In order to compare the numerical and analytical calculations, the results presented in Fig. 
2 are calculated for T = 0, A 31 = A 32 = 0, and vapor temperature T = 10 _3 i^ (i.e., in 
fact, for v z — 0). However, the atomic data (wavelengths, dipole moments and decay rates) 
are real. In particular, the spontaneous relaxation rates are not equal: 7 32 = 1. 66731. The 
intensities at the entrance to the medium are chosen so that u 2 = 0.59 • 10~ 4 . The dotted 
line in Fig. 2(b) is a calculation with formula $27\) for u\ = 0.55- 10~ 4 . One can see that the 
correspondence not only in the shape but also in quantitative characteristics of the process is 
excellent. The difference in parameter u\ Q is obviously due to the difference of spontaneous 
relaxation rates. As expected from analytical calculations, the transfer of energy to the 
THz radiation occurs with unity photon efficiency (the maximum value of It is exactly 
equal to 0.0021i3i(£ = 0)) at the length r max = 7.6 • 10 5 (( ma x = 5.87, this value is quite 
close to that calculated from Eq. (|26"D for u 2 ,^ = 0.59 • 10~ 4 : ( max = 6.25). Such a perfect 
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frequency conversion is due to the negligible decay of the dark state given by T and terms 
to the second order in / 7. For the parameters of Fig. 2, the maximum THz intensity 
It = O.OO2I/31 (£ = 0) corresponds to the Rabi frequency of only g? w 8 • 10 _4 73i. The 
behavior of the phase $ in Fig. 2(c) also follows the law cos$(C) = obtained analytically. 
The jumps in the phase occur at points where the intensity of the field being absorbed 
approaches zero, according to Eqs. (|| (b,d)). 

In reality, however, we do have both the Doppler broadening and the relaxation rate V 
of the coherence between states |1) and |2). Quantitatively, this considerably modifies the 
process. In Fig. 3 the spatial dependence of the field intensities and the phase $ are plotted 
for the vapor temperature T = 800 K (this gives the saturated vapor density of iV = 1.7 • 
10 15 cm -3 and corresponds to the most probable velocity of atoms of v p = 7.5 • 10 4 cm/ sec), 
T = 10 _4 73i (0.6 kHz), input Rabi frequencies g 3 i (r = 0) = 6O731, g 32 (t = 0) = 2O731 
(corresponding to laser intensities of J31 = 18.8 W/cm 2 and J32 = 2.1 W/ cm 2 ). We see that 
the dynamics of the intensities and the phase does not change qualitatively as compared 
to the case of negligible decay of the dark state (Fig. 2). The THz radiation is generated 
and reaches its maximum, It = 39 mW/cm 2 (Rabi frequency gT ~ 4.6 • 10 _3 73i) at the 
length r = 2.2 • 10 6 (for T = 800 K this corresponds to the real length of the gas cell of 
z = 4 cm) . A minor quantitative difference compared to the Fig. 2 case is that the optical 
length scale of the oscillations increases. This is simply a consequence of the optical length 
dependence on the total intensity, Eq. (p0|). However, the Doppler broadening, together with 
the rate T, considerably influence the efficiency of the frequency conversion. The efficiency 
reaches a maximum when the atoms are prepared in an almost dark state so that the energy 
dissipation is very weak. The dark state preparation relies in present case on the optical 
pumping. Therefore, the optical pumping rate should be much larger than the dark state 
decay in order to allow the population trapping in \NC). This requirement results in the 
following condition^: 

gl r 
— — — s> — 

7 2 + A 2 7' 

12 



where detuning A includes the Doppler shift: A = A31 — k^\v z pa A32 — ^32^- For the 
resonance A 31 = A 32 = and large Doppler broadening k 3 iv p ^> 7, this condition reduces 
to 

9o » - (k 3 iv P ) 2 (28) 

7 

or, in terms of intensity, 

x >> £ fk^V 16^70 
7 V 7 / 3A^ 

For the here considered case of 24 Mg atoms at temperature 800i^, the total input intensity 
should be larger than 3.75 ■ 10 3 ■ r/7 [W/cm 2 ]. An importance of the condition ( 28|) is 



demonstrated in Fig. 4. While the efficiency is approaching maximum: 77 = 2.06 • 10~ 3 , for 
the parameters of Fig. 3 where r/7 = 10~ 4 , it drops to 77 = 1.66 • 10~ 3 for r/7 = 2 • 10~ 3 in 
Fig. 4 (a). Moreover, the intensity oscillations decay quite fast, and at the length r > 4 - 10 6 
(C > 4.87) the optical fields are completely absorbed by the medium (not shown). Increase of 
the input intensity in Fig. 4 (b) leads to better trapping, weaker e.m. energy dissipation and, 
correspondingly, to larger conversion efficiency increasing to r] — 2.00 • 10~ 3 . Of course, with 
increasing input intensity the length scale of the conversion increases as well. Parameters of 
Fig. 4 (b) correspond to input intensities J31 = 470 W/cm 2 and J 32 = 52 W/cm 2 , and to the 
maximum THz intensity of It = 0.94 IV/ cm 2 (Rabi frequency gx ~ 2.3 • 10~ 2 7 3 i) generated 
at t = 1.1 ■ 10 7 (z = 20cm). 

Thus, a practical realization of the proposed here generation scheme requires large optical 
length (large atom density), fairly large input intensities of the optical fields and small 
decay rate T of the atomic coherence G\i. The rate V is determined by uncorrelated laser 
fluctuations, atomic collisions and other random phase disturbing processes. It has recently 
been experimentally demonstrated that, in a buffer gas, the coherence a 12 can survive a very 
large number of atomic collisions@@. For instance, the rate T < 50 Hz has been observed in 
experiment^. The correlation of the fluctuations in U31 and u; 3 2 radiations may be achieved 
if, for example, the 0^32 frequency is generated from cj 31 by operation of a Raman laser, as 
in Ref.0. 

13 



An important point in nonlinear frequency conversion is a phase matching. In our scheme, 
however, the optical frequencies are quite close in magnitude to each other. Therefore, 
the wave vectors k^ m change approximately equally due to dispersion of the medium, in- 
duced by the interaction with far-detuned states other than those of the A system. Thus, 
this phase mismatch may be neglected. Moreover, as it has been shown theoretically and 
experimentally^, the phase mismatch can be compensated for by very small value of two- 
photon detuning (A 3 i — A 32 ), if it does appear. 

In the particular case of Mg atoms, considered here as an illustrative example, the 
metastable states 3 3 Pi and 3 3 P2 are not the ground states. Therefore, one has to prepare 
an atomic ensemble in metastable states. This problem may be solved, for example, through 
a low-energy high-current electric discharge!. Another problem is that the state 3 3 Pi decays 
spontaneously to the ground state 3 1 S . However, the rate of this decay, 4.3 • 10 2 sec -1 = 
1.2 • 10 _5 73i is too slow to produce considerable effect on the process, presented here. Of 
course, all these problems may be avoided if one uses atoms with the metastable states |1) 
and 1 2) being the ground states, like in Sm atoms. 

In conclusion, we have proposed a scheme for generation of continuous-wave THz radia- 
tion. Our scheme is based on a nonlinear interaction of three e.m. waves in atomic media 
with a A configuration of levels. An important ingredient of the scheme is a preparation of 
atoms in the grey state. As a consequence, the e.m. fields propagate through the medium 
with very weak dissipation of energy. Therefore the process may be considered 
metric difference-frequency generation, where the THz frequency is the difference between 
two optical ones. We have obtained an analytical solution of the e.m. propagation problem 
which demonstrates that the intensities oscillate along the propagation path in a form of 
Jacobi elliptic functions. This allows us to predict that the photon conversion efficiency 
approaches unity in this technique, and to estimate the optical length at which the energy 
transfer from the optical field into the THz one is maximum. The analytical solution is 
confirmed by numerical calculations taking into account Doppler broadening and relaxation 
of the grey state. These calculations show that the efficiency of the THz generation remains 
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still very high in real situations, if the input optical intensity is sufficiently high to satisfy the 
condition (|29|). We note finally that the tuning of obtained THz radiation may be realized 
by shifting the metastable state energies via strong static magnetic or/and electric fields. 
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Figure captions 

Fig. 1. Closed A system with two metastable states |1) and |2). cu 31 and CU32 are the 
optical frequencies, ojt is the THz-range frequency. 

Fig. 2. Spatial variations of the optical (a) and THz (b) field intensities and the relative 
phase $ (c) in a vapor of 24 Mg atoms interacting with radiation in a closed A configuration 
of levels 3 3 Pi — 3 3 P2 — 4 3 S*i. For this system, the relaxation rates are 731 = 3.46 ■ 10 7 sec -1 , 
732 = I.66731, 721 = 2.6-10~ 14 73i, the wavelengths A31 = 517.27 nm, A 3i = 518.36 nm. Other 
parameters are: vapor temperature T = 10~ 3 K, T — 0, detunings A31 = A32 = 0, Rabi 
frequencies of input fields g 3 i (r = 0) = IO731, g 32 (r = 0) = O.I731 and gi 2 (r = 0) = 0. The 
dotted curve in (b) is a calculation with formula (|27|) for ul = 0.55 • 1CT 4 . 

Fig. 3. Spatial variations of the optical (a) and THz (b) field intensities and the relative 
phase $ (c) in a vapor of 24 Mg atoms for vapor temperature T = 800 K, T = 10~ 4 73i, 
detunings A31 = A 32 = 0, Rabi frequencies of input fields g%i (r = 0) = 6O731, g^ (r = 0) = 
2O731 and #i2 (t = 0) = 0. Other parameters are the same as in Fig. 2. 

Fig. 4. Spatial variations of the THz radiation intensity in a vapor of 24 Mg atoms for 
T = 2 • 10~ 3 73i and Rabi frequencies of input fields: (a) g%i (r = 0) = 6O731, #32 (r = 0) = 
2O731 and (b) g 3 i (r = 0) = 3OO731, (732 (r = 0) = IOO731. Other parameters are the same as 
in Fig. 3. Notice the different length scales in (a) and (b). 
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